RF electromagnetic field and vortex penetration in multilayered superconductors 
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Multilayered structure with a single superconductor layer and a single insulator layer deposited 
on a bulk superconductor is studied. General formulae for the vortex penetration-field of the su- 
perconductor layer and the magnetic field on the bulk superconductor which is shielded by the 
superconductor and insulator layers are derived with a rigorous calculation of the magnetic field 
attenuation in the multilayered structure. The formulae depend not only on the material and the 
thickness of the superconductor layer but also on the thickness of the insulator layer. The results 
can be applied to superconducting accelerating cavities with the multilayered structure. Using the 
formulae, a combination of the thicknesses of superconductor and insulator layers that can realize 
the enhanced rf breakdown field can be found for any given materials. 



Technologies to fabricate the superconducting rf cav- 
ities made of Nb have been advanced. The maximum 
accelerating gradient £^acc of the TESLA type 1.3 GHz 9- 
cell cavities during performance tests in vertical cryostats 
regularly exceed 35 MV/m at several laboratories. The 
gradient record had been increasing and recently two 9- 
cell cavities made from large grain Nb reached 45 MV/m 
at DESY [l|]. Further high gradients, however, would 
not be expected because their gradients are thought to 
be close to the empirical limit imposed by the thermo- 
dynamic critical field ~ 200 niT of Nb l^- A. Gurevich 
suggested [S, |j| that a multilayered nanoscale coating on 
Nb cavity may push up the rf breakdown field to the level 
of the vortex penetration-field of the coating materials at 
which the Bean-Livingston surface barrier p] disappears. 
While some experimental studies have been conducted on 
the subject based on the idea [6|,lZ|, not much theoretical 
progress followed on it. In fact, the best parameter set for 
the multilayer coating model such as thicknesses of lay- 
ers and choices of materials is not clear in the theoretical 
view. In this letter, the multilayered structure is carefully 
evaluated with a rigorous calculation on the electromag- 
netic field ditribution to keep its self-consistency. The 
resultant vortex penetration-field, the best combination 
of parameters, and materials are described. 

The multilayer coating model [3| consists of alternating 
layers of superconductor layers iS) and insulator layers 
{I). The simplest configuration with a single supercon- 
ductor layer and a single insulator layer is seen in Fig. [1] 
Each S layer is expected to withstand higher field than 
bulk Nb, and to shield the bulk Nb from the applied 
rf surface field Bq, because Bi (an rf surface field on 
the bulk Nb) is smaller than Bq. Then the multilayered 
structure is thought to withstand a higher field than the 
bulk Nb if Bq is smaller than the vortex penetration- 
fields of the top S layer and Bi is smaller than that of 
the bulk Nb. The vortex penetration-field of the S layer 
was given by B^ = 0o/47rA^ in the original paper [3J, 
where 0o = 2.07 x 10^^^ Wb is the flux quantum [1|, and 



A and ^ are a London penetration depth and a coher- 
ence length of the material of the S layer, respectively. 
This expression, however, has the same form as the vor- 
tex penetration-field of the semi-infinite superconductor, 
and does not depend on any parameters on the configu- 
ration of the multilayered structure such as the S layer 
thickness or the I layer thickness. In order to incorporate 
effects from the configuration of multilayered structure, 
we carried out rigorous caluculation on the distribution 
of magnetic field and Meissner current in the S layer. 

In order to derive the electromagnetic field in the 
multilayered structure, the Maxwell equations and the 
London equations should be solved with appropriate 
boundary conditions simultaneously. Contributions to 
the electromagnetic field distribution from the normal 
(unpaired) electrons of the superconductor and dielectric 
losses in the insulator are neglected. For simplicity, let 
us consider a model with a single S layer and a single 
T layer deposited on a bulk superconductor as shown in 
Fig. [TJ Table. U shows the parameters for the model, di 
is assumed to be zero or larger than a few nm to sup- 
press the Josephson coupling [J. All layers are parallel 
to the y-z plane and then perpendicular to the x-axis. 
The applied electric and magnetic field are assumed to 
be parallel to the layers. Further we assume the materials 
used for the S layer is extreme Type II superconductor 
Ai ^ ^1, and the S layer thickness is larger than the the 
coherence length dg >• ^i . Note that the S layer of our 
model is not necessarily a thin film hence the discussion 
below can be applied to any S layer with arbitrary thick- 
ness ds ^ Ci- Solving the Maxwell equations in the X 
layers, and the Maxwell-London equations in the S layers 
and in the bulk superconductor, we find 

Ai cosh ^H^ + (Aa + rfi) sinh ^H^ 
Bi = Bo— /:,. .t' , (1) 
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TABLE I. Parameters of the multilayered structure with a single superconductor layer and a single insulator layer deposited 
on a bulk superconductor. 



Region 



Material type 



Parameter 



I 

II 
III 



Superconductor layer 

Insulator layer 
Bulk superconductor 



Coherence length: ^i, London penetration depth: Ai (S> ^i), Thickness: ds (:§> ^i) 

Relative permittivity: tr, Thickness: dx (zero or larger than a few nm) 

London penetration depth: A2 




X (nm) 

FIG. 1. Examples of the magnetic field attenuations in the 
multilayered structure. Black lines show our formulae given 
above, and gray lines show the naive estimates with B — 
Boe~^''^^. Solid lines, dashed lines and dashed-dotted lines 
correspond to the material of the <S layer: NbN (Ai — ANbN = 
250 nm), MgBj (Ai = AmsB^ = 110 nm) and NbsSn (Ai = 
ANbsSn ~ 60 nm), respectively. The bulk superconductor is 
assumed to be Nb (A2 = ANb ~ 50 nm). The thickness of 
the S layer and the I layer are fixed at ds = 50 nm and 
dx — 20 nm. 



where i?i, Bu, and i?iii are magnetic fields in region I, 
II, and III, respectively. Note here that these equations 
are valid for dx ^ (i/e7A:)~^ ~ 1cm, vi^here k — oj/c, 
uj = 2Tr X 1.3 X 10^ s""'^ is the angular frequency and 
c is the speed of light. Since we consider an insulator 
thickness dx <S^ 1cm below, Eq. ([T]), ^ and (jS]) are ac- 
curate enough. It is easy to confirm that these equations 
are reduced to the well known expression for the semi- 
infinite superconductor given by B — Boe~^'^^ when 
the S layer and the bulk superconductor are the same 
material (Ai = A2) and the I layer vanishes {dx — > 0). 
Fig. [T] shows examples how a magnetic field attenuates 
in a multilayered structure. 

The vortex penetration-field can be evaluated by com- 
puting two forces acting on a vortex at a top of the 
S layer: a force from an image current of an image 
vortex which is introduced to satisfy a boundary con- 
dition of zero current normal to the surface, and an- 
other from a Meissner current jm due to existence of 
external field which can be computed from Eq. ^ with 
JMy — —{l/no)dBi/dx. Then the vortex penetration 
field is given by 
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(A2 + dx) cosh ^ 
which depends on both the S layer thickness ds and the 



I layer thickness dx- Note here that Eq. Q is reduced to 
the well-known expression (jyo/AnXi^i for the semi-infinite 
S layer {ds -> 00). As is obvious from Eq. (|4]), the vortex 
penetration-field decreases as the I layer thickness dx in- 
creases. This behavior can be understood from the above 
results that the magnetic field attenuates more rapidly in 
an S layer on a thick I layer than that on a thin X layer. 
This means that a Meissner current, which is propor- 
tional to a gradient of the magentic field, becomes larger 
in an S layer on a thick X layer than that on a thin I 
layer. As a result the vortex is strongly drawn into the 
S layer by the force from the Meissner current. 

Let us define the magnetic field attenuation ratio a by 
a = Bu/Bq. Then the shielded magnetic field Bu on the 
bulk Nb is given by aB^ when the applied magnetic field 
Bq equals B^. When the magnetic field attenuation in 
the S layer is enough for the shielded magnetic field aB,; 
to become smaller than 200 mT, which is thought to be 
the maximum field for the bulk Nb, the bulk Nb is safely 
protected. Then the maximum peak surface-field -Bmax 
is given by By. On the other hand when the magnetic 
field attenuation is not enough and aBy is larger than 
200 mT, Bmax is limited by a-^ x 200 mT. Thus we find 



-L^niax — 



By 



{aBy < 200 mT) 



a"ix200mT (aB„ > 200mT). 



(5) 



Fig. [2] shows contour plots of i?max- In order to improve 
Bmax, appropriate parameter regions should be chosen, 
otherwise -Bmax becomes smaller than that of the bulk 
Nb. As is obvious from comparison between two con- 
tour plots, the material of the S layer drastically changes 
the effect of the multilayered structure. As shown in 
Fig. [SJa), only small area of the parameter space is al- 
lowed for the enhancement of -Bmax with a single NbN 
layer. In addition, even if the best parameter set are cho- 
sen, the multilayered structure with a single NbN layer 
hardly can enhance -Bmax- On the other hand, as shown 
in Fig. [2ljb), NbsSn layer significantly increases -Bmax in 
large area of the parameter space because of its superior 
screening effect and its high vortex penetration-field. 

The general formulae of magnetic field distributions on 
structure of a couple of superconductor and insulator lay- 
ers on a bulk superconductor are solved self-consistently. 
Vortex penetration-field for a superconductor sheet with 
different field level on both the sides is also derived. Com- 
bining these results, we can survey optimum combination 
of thicknesses on S and I layers, and S layer material. 
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FIG. 2. Contour plots of the maximum peak surface-field 
Bmax- The abscissa represents the I layer thickness di and 
the ordinate represents the <S layer thickness ds- Values writ- 
ten in the plot area are Bmax in the unit of mT. The top 
and bottom figures correspond to materials of the S layer: 

(a) NbN (Ai = ANbN = 250 nm, ^i = ^NbN = 5nm) and 

(b) NbsSn (Ai = ANbaSn = 60 nm, ^i = ^NbgSn = 3nm), 
respectively. The bulk superconductor is assumed to be Nb 
(A2 = ANb — 50 nm). 



These procedure shows how the alternative layered struc- 
ture can contribute to raise the maximum rf surface field. 
Hence it is possible for us to evaluate how the introduc- 
tion of this structure on the inner-surface of supercon- 
ducting cavity improves its achievable accelerating gra- 
dient. 
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